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ABSTRACT 

We  discuss  a  reliability  model  which  reflects  the  dynamic  dependency  between 
system  failure  and  system  stress.  In  particular,  a  shot-noise  process  Is  used  to 
model  "residual  system  stress,  ",%vhlch  In  turn  drives  a  doubly  stochastic  Poisson 
process  model  for  system  failures.  Intuitively,  residual  stress  (or  susceptibility  to 
failure)  may  vary  In  a  random  manner  yet  be  essentially  unobservable,  while 
system  failures  may  be  readily  detectable  and  observable.  Shot-noise 
distributions  have  a  richness  and  subtlety  which  suggest  untapped  potential  for 
applications.  The  doubly  stochastic  Poisson  process  provides  a  reasonable 
framework  for  modeling  randomly  varying  rates  of  occurrence  In  a  broad  variety 
of  settings. 
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1.  INTRODUCTION 


Tills  note  Is  motivated  by  the  general  approach  to  failure  modeling  discussed 
in  Lemolne  and  Wenocur  [4].  Under  this  approach,  system  state  or  wear  and  tear 
Is  modeled  by  an  appropriately  chosen  random  process;  and  the  occurrence  of 
fatal  trauma  Is  modeled  by  a  Poisson  process  whose  rate  function  Is  state 
dependent.  The  system  Is  said  to  fall  when  either  wear  and  tear  accumulates 
beyond  an  acceptable  or  safe  level  or  fatal  trauma  occurs.  The  rate  at  which 
fatal  trauma  can  occur  depends  on  the  state  of  the  system.  The  approach  has 
Interesting  implications.  It  provides  revealing  Insights  Into  most  of  the  famous 
and  frequently  used  lifetime  distributions  In  reliability  theory.  Including  the 
Makeham,  Gompertz,  Welbull,  Rayleigh  and  Gumbel  distributions;  In  fact,  these 
classic  models  are  obtained  In  a  unified  and  straightforward  manner.  Moreover, 
the  approach  suggests  Intuitively  appealing  and  computationally  tractable  ways 
of  enhancing  these  standard  failure  models  and  for  developing  new  ones. 

In  [4]  detailed  examples  of  the  failure  modeling  approach  are  provided  using 
deterministic  models  and  diffusion  process  models  for  system  state;  use  of  a  shot- 
noise  model  for  system  state  Is  discussed  In  fairly  general  terms.  In  the  present 
note  we  provide  a  more  detailed  discussion  of  the  approach  based  on  using  a 
shot-noise  model  for  system  stress.  And  we  assume  that  the  rate  of  occurence  of 
fatal  trauma  Is  proportional  to  "residual  system  stress".  The  time  to  system 
failure  Is  then  the  epoch  of  the  first  count  In  a  (doubly  stochastic)  Poisson  process 
whose  rate  function  Is  the  shot-noise  process. 

In  using  a  shot-noise  model  for  system  stress  we  assume  that  the  system  Is 
subjected  to  "shots"  or  Jolts  according  to  a  stochastic  point  process.  A  Jolt  may 
consist  of  an  internal  component  malfunctioning  or  an  external  "blow"  to  the 
system.  Jolts  Induce  stress  on  the  system  when  they  occur.  However,  If  the 
system  survives  the  Jolt  it  may  recover  to  some  extent.  For  Instance,  a  sudden 
and  unexpected  surge  of  power  In  the  circuit  of  a  control  system  may  temporarily 
Increase  the  likelihood  of  system  failure,  but  the  overload  Itself  decays  rapidly. 
For  another  example,  the  mortality  rate  for  persons  who  have  suffered  a  heart 
attack  declines  with  the  elapsed  time  since  the  trauma.  In  this  case,  the  heart 
actually  repairs  Itself  to  a  degree. 

The  remainder  of  this  note  Is  organized  as  follows.  In  the  next  section. 
Section  2,  we  derive  a  class  of  failure  distributions  using  the  shot-noise  model  for 
system  stress,  and  discuss  properties  and  structure.  Following  this,  we  discuss  In 
Section  3  the  shot-noise  process  and  a  related  stochastic  Integral  representation; 
this  parallels  the  technical  development  In  [4)  based  on  stochastic  differential 
equations  and  diffusion  processes. 

For  background  on  shot-noise  processes  and  distributions  see  Rice  [5],  Cox 
and  Isham  (2),  and  Bondesson  [lj.  Doubly  stochastic  Poisson  processes  are 
discussed  by  Grandell  [3].  The  gamma  process  Is  discussed  In  Takacs  [6). 
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2.  SHOT-NOISE  AND  FAILURE  DISTRIBUTIONS 

Suppose  that  a  system  Is  subjected  to  shots  or  Jolts  according  to  a  Poisson 
process  with  deterministic  rate  function  {m(u)  ,  u^O}.  Suppose  that  If  a  Jolt  of 
magnitude  D  occurs  at  epoch  S  then  at  time  S+t  the  contribution  of  the  Jolt 
to  the  system  stress  Is  Dh(t),  where  h  Is  a  non-negative  function  which  tends 
to  zero  as  t  -*°°  .  (  The  function  h  Is  also  assumed  to  vanish  on  (— 00 , 0)  and 
to  be  lntegrable  over  finite  Intervals.)  In  other  words,  shot-induced  stress  Is 
additive  and  decays  with  time  according  to  the  rate  function  h.  If  {Sn  ,  n^l} 
are  the  epochs  of  shot  occurrences  and  {D„  ,  n25l}  are  the  magnitudes  of  the 
successive  Jolts,  then  the  "residual  system  stress"  at  time  t,  say  X(t)  ,  Is  given  by 

X(t)  -  2D„h(t-S„)  (l) 

n  “  1 


where  we  take  X(0)=0  .  With  stress  comes  susceptibility  to  failure,  and  It  Is 
assumed  that  the  occurrence  rate  of  fatal  trauma  Is  an  Increasing  function  of 
residual  system  stress.  In  particular,  If  T  Is  the  time  to  system  failure,  we 
assume  that  P{  TS8+u  |  T>u  ,  X(u)  *=  x  }  =  8kx  +  o(8)  for  each  time 
point  u  and  some  parameter  k  ,  so  that  kx  can  be  Interpreted  as  the  rate  of 
occurrence  of  fatal  trauma  when  the  level  of  system  stress  Is  x.  Thus,  the  failure 
time  T  Is  the  epoch  of  the  first  count  In  a  Poisson  process  with  rate  function 
|kX(u)  ,  u^O},  and  therefore 

t 

P{  T  >  t  }  =  E[  exp[-JkX(u)du]  ].  (2) 

o 

The  usual  choice  for  the  attenuation  or  recovery  function  h  In  the  literature 
Is  the  function  h(u)=5cxp(—au)  If  u  2s  0  and  h(u)  =*  0  If  u  <  0  ,  where  a  Is 
some  positive  parameter,  that  Is,  exponential  decay.  Alternate  choices  for 
modeling  the  pattern  of  shot  occurrences  Include  a  renewal  process,  a  semi- 
Markov  process,  or  a  cluster  point  process.  In  this  note  we  assume  a  Poisson 
pattern  for  shot  occurrences.  We  also  take  the  shot  occurrences  and  Jolt 
magnitudes  to  be  Independent  of  one  another,  and  the  successive  Jolt  magnitudes 
to  be  Independent  and  distributed  as  a  random  variable  D  which  Is  not 
Identically  0  and  has  a  proper  distribution. 

We  now  compute  the  right  side  of  equation  (2)  and  related  quantities,  discuss 
properties  of  the  resulting  formulas,  and  then  look  at  some  examples. 

To  begin,  let  N(t)  be  the  number  of  shot  occurrences  up  to  time  t  , 

t 

M(t)  “  J*m(u)du  (3) 

o 


and 
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Y(t)  =  /kX(u)du. 
o 


We  assume  that  E{N(t)}  -  M(t)  Is  finite  for  each  t  ,  from  which  It  follows 
that  Y(t)  Is  finite  with  probability  one.  Put 


H(t)  =  jh(u)du. 


If  N(t)  =  n,  then 


Y(t)  =  2kDjH(t— Sj) 

i-i 


and  the  shot  epochs  Sj,  *  *  •  ,Sn  are  Independent  and  distributed  as  a  random 
variable  U  with  density  m(u)/M(t)  for  O^u^t.  Thus,  If  L  Is  the  Laplace 
transform  of  the  distribution  of  D  ,  then  for  s^O 

E{exp[-sY(t)]|N(t)=n}  =  (e{  exp[-skDH(t-U)]  })"  (7) 

=  |/L(skH(t-U))[m(uVM(t)]du|  . 

Unconditioning  on  the  value  of  N(t)  ,  we  find  that 

t 

E{  exp[-sY(t)]}  =  exp(-M(t)]exp[/L(skH(u))m(t-u)du].  (8) 

o 

Putting  3=1  In  equation  (8)  gives 


P{  T>t  }  =  exp[— M(t)]  exp[/L(kH(u))m(t—  u)duj.  (9) 

o 

If  M(t)  -* 00  as  t  -♦ 00  then  equation  (9)  defines  a  proper  probability 
distribution  concentrated  on  (0,»).  We  assume  henceforth  that  this  Is  the  case. 
Let  R(t)  denote  the  quantity  on  the  right  side  of  equation  (9).  If  the  function 
m(«)  driving  the  process  of  shot  occurrences  Is  almost  everywhere  differentiable, 
then  the  distribution  of  T  has  density  R(t)r(t)  and  failure  rate  r(t),  where 

t 

r(t)  =  m(t)  -  m(0)L(kH(t))  -  jL(kH(u))m'(t-u)du.  (10) 
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Now  consider  N(T),  the  number  of  Jolts  leading  to  system  failure.  We 
Indicate  briefly  how  to  obtain  the  generating  function  of  N(T).  We  begin  with 
the  Joint  distribution  of  T  and  N(T).  Consider  the  difference 

P{t-h<TSt,  N(t)=n}  =  P{T>t-h,N(t)=n}-P{T>t,N(t)=n}  (11) 

Formally,  If  we  divide  equation  (11)  by  h  and  let  h  tend  to  0,  we  obtain  the  Joint 
density  of  T  and  N(T).  This  limiting  process  can  be  rigorously  Justified  by 
appealing  to  the  absolute  continuity  of  the  distribution  of  T.  We  now  proceed  by 
expanding  the  left-hand-side  of  (11)  Into  more  tractable  quantities.  Standard 
probabilistic  arguments  lead  to: 

P{T>t-h  ,  N(t)=n}  -  P{T>t,N(t)=n}  =  g(t-h,n)-g(t,n)  (12) 

-g(t-h,n)[l-P{N(t)-N(t-h)=0}]  +  g(t-h,n-l)P{N(t)-N(t“h)=l}  +  o(h), 

where  g(t,n)  =  P{T>t  ,  N(t)=n}.  Dividing  equation  (12)  through  by  h  and 
letting  h-*0  gives 

P(Tcdt  ,  N(T)=n)  =  [g(t,n— 1)—  g(t,n)]m(t)dt  —  g'(t,n)dt.  (13) 

Now 


g(t,n)  =  [ J L(kH(u))m(t -u)du] "exp [  — M(t)]/n !  ,  (14) 

o 

and  substituting  (14)  Into  equation  (13)  yields 

P(T€dt  ,  N(T)=n)  =  g(t,n— l)r(t)dt .  (15) 

A  simple  Integration  and  summation  now  show  that  for  0  S  0  2S  l, 

«  t 

E[0n(t)]  =  j0exp[0/L(kH(u))m(t— u)du]exp[— M(t)]r(t)dt.  (16) 

o  o 

Finally,  differentiating  (18)  with  respect  to  0  and  setting  0  =  1  gives 

E[N(T)]  =  1  +  jf|/L(kH(u))m(t-u)du|R(t)r(t)dt.  (17) 


Other  quantities  of  interest  are  the  mean  and  variance  of  Y(t). 
Differentiating  the  tranform  (8)  with  respect  to  a  and  then  setting  a  =  0,  we 
see  that 

t 

E{Y(t)}  =  E{D}  /kH(u)m(t-u)du,  (18) 

o 

and  If  E{Y(t)}  is  finite, 

t 

Var{Y(t)}  =  E{D2}  J[kH(u)l2m(t-u)du.  (19) 

o 

The  most  tractable  circumstance  occurs  when  m(«)  Is  a  step  function  or 
constant.  In  which  case  equation  (10)  reduces  to 

r(t)  =  m(t)  —  m(0)L(kH(t)).  (20) 

Since  this  case  covers  a  multitude  of  possibilities,  we  assume  henceforth  that  the 
rate  function  driving  the  shot  process  Is  either  a  step  function  or  constant.  To 
preclude  tedious  nuisance,  we  assume  that  m(0)  *  m(0+)  >  0.  We  now 
consider  some  specific  examples. 

Example  1:  Suppose  that  the  Jolt  magnitudes  are  constant  and  that  the  Jolt- 
Induced  stress  decays  slowly.  In  particular,  suppose  that  D  ®  d  and  that 
h(u)=l/(l  +  u)  for  u  2:  0  .  Then 

r(t)  *  m(t)  -  [m(0)/(l  +  t)kd],  (21) 

R(t)  =  exp[-M(t)](l  +  t)m<°)  (22) 

If  kd  =  1  ,  and 

R(t)  =  exp[—  M(t)]exp[m(0){(l  +  t)1-kd  —  l}/(l  —  kd)]  (23) 
otherwise.  If  M(t)=mt,  then  (16)  reduces  to 

R(t)  =*  exp[— mt](l  +  t)m.  (24) 

Example  2:  Suppose  that  D  has  a  geometric  distribution;  specifically, 
suppose  that  P{  D  *  n  }  "  qpn  for  n  *  0,1,2,...  ,  where  0  <  p  <  1 
Assume  that  stress  decays  slowly  as  In  Example  1,  and  that  k  =  1  .  Then 
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r(t)  =  m(t)  -  m(0)q[(l  +  t)/(q  +  t)]  (25) 


and 


R(t)  =  exp[m(0)qt  -  M(t)]  {(q  +  t)/q}  m^°^qp.  (26) 

Example  3:  Suppose  there  Is  a  delayed  reaction  to  a  Jolt,  followed  by  slow 
decay;  In  particular,  assume  that  h(v)  =  ln(l  +  av)/(l  +  av)  for  v2:0,  where 
a  Is  some  positive  parameter.  Further,  suppose  that  the  distribution  of  D  has  a 
heavy  tall;  In  particular,  assume  the  distribution  of  D  has  transform 
exp[“’b(2s)1/2]  for  s  2s  0  where  b  is  some  positive  parameter.  Putting 
c  —  ba_  112  and  k  =  1  gives 

r(t)  =  m(t)—  m(0)  [1/(1  +  at)6],  (27) 

which  Is  very  similar  to  (21)  .  If  also  M(t)  —  mt  and  c  =  1  then 

R(t)  =  exp[— mt](l  +  at)m.  (28) 


Example  4:  Suppose  that  shot-induced  stress  decays  at  an  exponential  rate;  In 
particular,  assume  h(v)  *  exp[— av]  for  vS:0,  where  a  Is  some  positive 
parameter.  Further,  suppose  that  D  has  an  exponential  distribution  with 
parameter  b.  Then 

r(t)  —  m(t)  —  m(0)  (ab/(ab  +  k  —  kexp[— at])}  (29) 

and 

R(t)  =  exp  [{m(0)abt/(ab  +  k)}  —  M(t)]« 

{(ab  +  k  -  kexp[— at])  /  ab}  rn(°)b/C*b+k)  (30) 

In  particular,  If  M(t)  =  mt  and  k  =  1,  then 

R(t)  =  exp[-mabt/(l  +  ab)]  {(1  +  ab  -  exp[“at]yab}mb/^  +  ab).  (31) 

In  this  last  case,  residual  system  stress  X(t)  has  a  limit  distribution  (as  t  -*  “) 
which  Is  gamma  with  location  parameter  b  and  shape  parameter  m/a  ;  cf. 
Bondesson  [l].  This  is  a  remarkable  result  In  that  a  gamma  distribution  with 
arbitrary  shape  parameter  arises  In  a  physical  model. 


S.  A  GAMMA  PROCESS  DRIVING  SYSTEM  STRESS 

Let  {Z(u)  ,  u^O}  be  a  compound  Poisson  process  with  Jump-rate  function 
(m(u)  ,  u&0}  and  Jump  sizes  distributed  as  D.  Note  that  the  series 
representation  for  X(t)  given  by  (1)  truncates  at  N(t),  the  number  of  Jolts  up  to 
time  t.  Thus,  X(t)  admits  the  equivalent  representation 


x(<0  *  fh(t-u)Z(du)  (32) 

o 

where  the  right  side  of  (32)  Is  a  Rlemann-Stleltjes  Integral.  This  suggests 
replacing  the  compound  Poisson  process  Z(*)  to  drive  system  stress  by  some 
other  process  which  also  has  Independent  and  positive  Increments.  In  particular 
let  {Z*(u)  ,  uStO}  be  a  gamma  process  (cf.  (6))  for  which 

E{exp[— sZ*(u)]}  =  {h/(b+s)}u  (33) 

for  all  u  ,  s  2  0  ,  where  b  Is  some  positive  parameter.  (The  Increments  of  Z*(.) 
are  in  fact  stationary.)  Now  let  4 

t 

X*(t)  =  /h(t-u)Z*(du)  (34) 

0 

and  assume  that  the  residual  system  stress  at  time  t  Is  expressed  by  the 
Rlemann-Stleltjes  Integral  on  the  right  side  of  (34).  Further,  suppose  that  the 
system  failure  time  T  Is  the  epoch^  of  the  first  count  in  a  Poisson  process  with 
rate  function  {kX  (u)  ,  uS:0}.  If  R*(t)  =  P{T*>t}  then 

t 

R*(t)  “  E{exp[-/kX*(u)du]}  (35) 

0 

To  compute  the  right  side  of  (35),  we  first  note  that 

t  t 

J*X*(u)du  =  /H(t-u)Z’(du),  (36) 

0  0 

where  H(.)  is  given  by  (5).  ir  J  Is  a  positive  Integer  then 

ln(E{exp{— k2H(t-(jt/J))[Z*Qt/J)— Z*((j-l)t/J)]}})  =  (37) 

2  (t/J)ln{b/[b + kH(t — (j  t/J)]}. 
j-i 
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